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Abstract 

We review the behavior of epidemic spreading on complex networks 
in which there are explicit correlations among the degrees of connected 
vertices. 

1 Introduction 

Complex networks arising in the modeling of many social, natural, and tech- 
nological systems are often growing and self-organizing objects characterized 
by peculiar topological properties ^ Hj. Many empirical evidences have 
prompted that most of the times the resulting network's topology exhibits 
emergent phenomena which cannot be explained by merely extrapolating 
the local properties of their constituents. Among these phenomena, two 
of them appear ubiquitous in growing networks. The first one concerns the 
small-world property [HI, that is defined by an average path length — average 
distance between any pair of vertices — increasing very slowly (usually loga- 
rithmically) with the network size, A^. The second one finds its manifestation 
in the scale-free (SF) degree distribution This implies that the probabil- 
ity P{k) that a vertex has degree k — it is connected to k other vertices — is 
characterized by a power-law behavior P{k) ~ k^"' , where 2 < 7 < 3 is a 
characteristic exponent. 

The statistical physics approach ^ [2] has been proved a very valuable 
tool for the understanding and modeling of these emergent phenomena in 



growing networks and has stimulated a more detailed topological characteri- 
zation of several social and technological networks. In particular, it has been 
recognized that many of these networks possess non-trivial degree correla- 
tions [31 El d . The use of statistical physics tools has also evidentiated 
several surprising results concerning dynamical processes taking place on 
top of complex networks. In particular, the absence of the percolation [SI IHl 
and epidemic ^1 ^1 El E! thresholds in uncorrelated scale-free 
(SF) networks has hit the community because of its potential practical im- 
plications. The absence of the percolative threshold, indeed, prompts to an 
exceptional tolerance to random damages JH]- On the other hand, the lack 
of any epidemic threshold makes SF networks the ideal media for the propa- 
gation of infections, bugs, or unsolicited information jj^. While the study of 
uncorrelated complex networks is a fundamental step in the understanding 
of the physical properties of many systems T "71, yet correlations may dras- 
tically change the obtained results, as several recent works addressing the 
effect of correlations in epidemic spreading have shown ^1 El 1201 [22] . 

Here we want to provide a review of recent results concerning the epi- 
demic spreading in random correlated complex networks. We will consider 
the susceptible-infected-susceptible (SIS) and susceptible-infected-removed 
(SIR) models |23l 1241 125j and we will provide an analytical description that 
includes two vertices degree correlations in the dynamical evolution of the 
infection prevalence. This will allow us to relate the presence or absence of 
epidemic threshold to the eigenvalue spectra of certain connectivity matrices 
of the networks. In particular, in the case of scale-free networks it is possible 
to show that for the SIS model, a SF degree distribution P{k) ~ with 
2 < 7 < 3 in unstructured networks with any kind of degree correlations 
is a sufficient condition for a null epidemic threshold in the thermodynamic 
limit. For the SIR model, the same sufficient condition applies if the mini- 
mum possible degree of the graph is kmin > 2. The SIR model with kmm = 1 
has always a null threshold unless the SF behavior is originated only by min- 
imum degree vertices. In other words, under very general conditions, the 
presence of two-point degree correlations does not alter the extreme weak- 
ness of SF networks to epidemic diffusion. The present results are derived 
from the divergence of the nearest neighbors average degree , which stems 
from the degree detailed balance condition , to be satisfied in all physical 
networks. 

2 Correlated complex networks 

In the following we shall consider unstructured undirected networks, in 
which all vertices within a given degree class can be considered statistically 
equivalent. Thus our results will not apply to structured networks in which a 
distance or time ordering can be defined; for instance, when the small-world 
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property is not present [26| |27| . We will consider in particular the subset of 
undirected Markovian random networks |5J , that are completely defined by 
the degree distribution P{k) ant the conditional probability P(k' \ k) that a 
vertex of degree k is connected to a vertex of degree k' . These two functions 
can have any form and are assumed to be normalized, i.e. 

Y^p{k) = Y,p{k'\k) = i. (1) 

k k' 

The term "Markovian" refers to the fact that, in our approximation, all 
higher-order correlation functions can be obtained as a combination of the 
two fundamental functions P{k) and P(k' \k). In fact, this approximation 
represents a natural step towards a more complex description of real net- 
works. In this sense, the Erdos-Renyi (ER) model j2H] (defined starting 
from a set of N vertices that are connected in pairs with an independent 
probability p) can be viewed as the zero-th order approximation, where the 
average degree is the only fixed parameter. The ER model is thus defined 
by the ensemble of all possible networks with a given average degree, but 
completely random at all other respects. The first-order approximation has 
been recently introduced by realizing that many real-world networked sys- 
tems possess a more complex degree distribution than that predicted by the 
ER model (a Poisson distribution |22|). In this approximation the whole 
degree distribution, P{k), is chosen as a constrain whereas the rest of prop- 
erties are left at random |3n| 13 H I32[ I17j . Even though this approximation 
represents a quantitative step forward, it only takes into account local prop- 
erties and, therefore, it neglects possible correlations among different ver- 
tices, correlations that, on the other hand, are present in real networks j^[7j. 
Thus, it is quite natural to introduce the second-order approximation as that 
with a fixed degree distribution, P{k), and a two-point correlation function, 
P{k' I k), but totally random to all other respects. As we will see, in this 
case the approximation must be carefully made since, due to the two-point 
correlation constrain, the fundamental functions P{k) and P{k' \ k) must 
satisfy a peculiar detailed balance condition. 

The degree distribution usually identifies two kinds of networks. A first 
class, which includes classical models of random graphs [2Hj, is characterized 
by an exponentially bounded degree distribution. A second one refers to SF 
networks in which the degree distribution takes the form P{k) ~ Ck~"', 
usually with 2<7<3^|5]. In this case the network shows a very high 
level of degree heterogeneity, signalled by unbounded degree fluctuations. 
Indeed, the second moment of the degree distribution, (A:^), diverges in 
the thermodynamic limit — > oo, where kc is the maximum degree of the 
network. It is worth recalling that, in growing networks, kc is related to 
the network size N as kc ^ /v^i/(7-i) j^J. Noticeably, it is the large degree 
heterogeneity of SF networks that is at the origin of their extreme weakness 
towards epidemic spreading. 
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2.1 Assortative and disassortative mixing 



A direct study of the conditional probability P{k' \ k) in data from real 
networks usually yields results that are very noisy and difficult to interpret. 
In order to characterize the degree correlations, it is more useful to work 
with the average nearest neighbors degree (ANND) of the vertices of degree 
k 0], defined by 



k' 

and to plot it as a function of the degree k. When two-point correlations 
are not present in the network, the conditional probability takes the form 



P^'={k' I k) = k'P{k')/ (k), and the ANND reads k^^^ik) = (fc^) / (A;), which 



is independent on k. On the contrary, an explicit dependence of knn{k) on k 
necessary implies the existence of non-trivial correlations, as often measured 
in real networks [UEj. For instance, in many social networks it is observed 
that vertices with high degree connect more preferably to highly connected 
vertices; a property referred to as "assortative mixing". In this case, knn{k) 
is an increasing function of k. On the opposite side, many technological and 
biological networks show "disassortative mixing" ; i.e. highly connected ver- 
tices are preferably connected to vertices with low degree and, consequently, 
knn{k) is a decreasing function of k. Then, the ANND provides an easy and 
powerful way to quantify two-point degree correlations, avoiding the fine 
details contained in the full conditional probability P{k' \k). 

2.2 Degree detailed balance condition 

A key relation holding for all physical networks is that all edges must point 
from one vertex to another. This rather obvious observation turns out to 
have important implications since it forces the fundamental functions P{k) 
and P{k' \ k) to satisfy the following degree detailed balance condition |25 



This condition states that the total number of edges pointing from vertices 
with degree k to vertices of degree k' must be equal to the total number of 
edges that point from vertices with degree k' to vertices of degree k. This 
relation is extremely important since it constraints the possible form of the 
conditional probability P{k' \ k) once P{k) is given. It may be surprising 
that such a detailed balance condition exists since, in fact, networks are the 
result of a multiplicative random process and, in principle, detailed balance 
conditions only holds for systems driven by additive noise [33]. In fact, the 
usual detailed balance condition is the same as Eq. Q without the prefactors 
k and k' . These very prefactors account for the multiplicative character of 
the network and Eq. (jSJ can be viewed as a closure condition that guarantees 
the existence of the network. There is a simple way to derive this condition. 




(2) 




(3) 
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Let Nk be the number of vertices with degree k. Obviously, Nk = N 
and, consequently, we can define the degree distribution as 



P{k) = f . (4) 

The function P{k) alone does not define completely the topology of the 
network, because it says nothing about how vertices are connected to each 
other. Thus, we need to define additionally the matrix of connections among 
vertices of different degrees. Let iV^ be a symmetric matrix measuring the 
total number of edges between vertices of degree k and vertices of degree k' , 
when k ^ k\ and two times the number of self-connections, when k = k' . It 
is not difficult to realize that this matrix fulfills the identities 



k' 



= kNk. (5) 
EE^'^.'^' = WA^- (6) 

k k' 

The first of this relations simply states that the number of edges emanating 
from all vertices of degree k is kNk , while the second indicates that the sum 
of all the vertices's degrees is equal to two times the number of edges. The 
identity © allows us to define the joint probability 

P(M') = |^. (7) 

where the symmetric function (2 — 5k,k')P{k, k') is the probability that a 
randomly chosen edge connects two vertices of degrees k and k' . The cor- 
relation coefficient computed from this joint probability has been recently 
used in Ref. |7| in order to quantify two-point degree correlations. The tran- 
sition probability P{k' \k), defined as the probability that an edge from a k 
vertex points to a k' vertex, can be easily written as 

from where the detailed balance condition arises as a consequence of the 
symmetry of P{k, k') (or Nk^k')- 

From the degree detailed balance condition it is possible to derive some 
general exact results concerning the behavior of A^nn(^! ^c) and of (knri) u — 

P{k)knn{k, kc) in SF networks |22I- In these two functions we have now 
made explicit the kc dependence originated by the upper cut-off of the k 
sum and that must be taken into account since it is a possible source of 
divergences in the thermodynamic limit. The results that we will derive will 
turn out to be fundamental in determining the epidemic spreading behavior 
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in these networks. Let us start by multiplying by a A; factor both terms of 
Eq. Q and summing over k' and k. We obtain 



(^k') = Y,k'P{k')Y,kP{k\k% (9) 

k' k 

In SF networks with 2 < 7 < 3 we have that the second moment of the 
degree distribution diverges as {k"^) r^ki ^ *. We thus obtain from Eq. ©, 
using the definition (|2)). 

V k'Pik')~knnik', kc) ^ (10) 

where C is the constant prefactor from the degree distribution. In the case 
of disassortative mixing [7], the function knn{k',kc) is decreasing with k' 
and, since k'P{k') is an integrable function, the l.h.s. of Eq. ()10() has no 
divergence related to the sum over k' . This imphes that the divergence must 
be contained in the kc dependence of knn{k' , kc)- In other words, the function 
knnik',kc) — > 00 for A;c ^ cxD in a non-zero measure set. In the case of 
assortative mixing, knnik' , k^) is an increasing function of k' and, depending 
on its rate of growth, there may be singularities associated to the sum over 
k' . Therefore, this case has to be analyzed in detail. Let us assume that the 
ANND grows as knn{k', kc) ~ ak'f^, /? > 0, when k' ^00. If /? < 7 — 2, again 
there is no singularity related to the sum over k' and the previous argument 
for disassortative mixing holds. When 7 — 2 < /? < 1 there is a singularity 
coming from the sum over k' of the type a/cf "^K However, since Eq. (jlUj) 
comes from an identity, the singularity on the l.h.s. must match both the 
exponent of kc and the prefactor on the r.h.s. In the case 7 — 2 < /3 < 1, the 
singularity coming from the sum is not strong enough to match the r.h.s. 
of Eq. (fTU]) since /3 — (7 — 2) < 3 — 7. Thus, the function knnik' , kc) must 
also diverge when fee ^ 00 in a non-zero measure set. Finally, when (3 > 1 
the singularity associated to the sum is too strong, forcing the prefactor to 
scale as a ~ rkl ^ and the ANND as knnik', kc) — rkl ^k'^. It is easy 
to realize that r < 1, since the ANND cannot be larger than kc- Plugging 
the knnik' , kc) dependence into Eq. ([T?H) and simplifying common factors, 
we obtain the identity at the level of prefactors 

r 1 

(11) 



2-7 + /3 3-7 

Since /3 > 1 and r < 1, the prefactor in the l.h.s. of Eq. (|1H) is smaller than 
the one of the r.h.s. This fact implies that the tail of the distribution in the 
l.h.s. of Eq. (|lfl|) cannot account for the whole divergence of its r.h.s. This 

*For 7 = 3 the second moment diverges as (fc^) ^ Inkc but the argument, though more 
involved, is still valid. 
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means that the sum is not the only source of divergences and, therefore, the 
ANND must diverge at some other point. 

In summary, the function knnik',kc) must diverge when fee — > oo in a 
non-zero measure set independently of the correlation behavior. The large 
kc singularity of the ANND can then be used to evaluate the quantity 

(J^nri) ^ — P{k)knnik, kc) , (12) 
k 

where we have explicitly considered kc as a growing function of the network 
size N. The r.h.s. of this equation is a sum of positive terms and diverges 
with kc at least as knn{k, kc) both in the disassortative or assortative cases. 
In other words, all SF networks with 2 < 7 < 3 must present a {knn) —>■ 
00 for N — > 00. This statement is independent of the structure of the 
correlations present in the network. The quantity (knn) is particularly 
useful in model analysis and real data measurements. Degree correlation 
functions can be measured in several networks, but measurements are always 
performed in the presence of a finite kc that does not allow to exploit the 
singularity of the function knn{k,kc). The most convenient way to exploit 
the infinite size singularity is therefore to measure the (knn) for increasing 
network sizes. 

The quantity (knn) is very important in defining the properties of 
spreading processes in networks since it measures the number of individuals 
that can be infected in a few contagions. We shall discuss this point in 
relation to some specific epidemic models in the next sections. 

3 The SIS model 

As a first prototypical example for examining the properties of epidemic dy- 
namics in SF networks we consider the susceptible-infected-susceptible (SIS) 
model j23j , in which each vertex represents an individual of the population 
and the edges represent the physical interactions among which the infection 
propagates. Each individual can be either in a susceptible or infected state. 
Susceptible individuals become infected with probability A if at least one 
of the neighbors is infected. Infected vertices, on the other hand, recover 
and become susceptible again with probability one. A different recovery 
probability can be considered by a proper rescaling of A and the time. This 
model is conceived for representing endemic infections which do not confer 
permanent immunity, allowing individuals to go through the stochastic cycle 
susceptible infected susceptible by contracting the infection over and 
over again. 
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3.1 Uncorrelated homogeneous networks 

In uncorrelated homogeneous networks, in which each vertex has more or 
less the same number of connections, A; ~ {k), a general result states the 
existence of a finite epidemic threshold, separating an infected (endemic) 
phase, with a finite average density of infected individuals, from a healthy 
phase, in which the infection dies out exponentially fast j25j. This is for 
instance the case of random networks with exponentially bounded degree 
distribution. 

This result can be recovered by considering the dynamical evolution of 
the average density of infected individuals p{t) (the prevalence) present in 
the network. The SIS model in homogeneous uncorrelated networks at a 
mean-field level is described by the following rate equation 



In this equation we have neglected higher order terms, since we are inter- 
ested in the onset of the endemic state, close to the point p{t) r-^ 0. Also, 
we have neglected correlations among vertices. That is, the probability of 
infection of a new vertex — the second term in Eq. H13|l — is proportional to 
the infection rate A, to the probability that a vertex is healthy, 1 — p{t), 
and to the probability that a edge in a healthy vertex points to an infected 
vertex. This last quantity, assuming the homogeneous mixing hypothesis^ 
is approximated for homogeneous networks as {k) p{t), i.e. proportional to 
the average number of connections and to the density of infected individu- 
als. From Eq. (|13j) it can be proved the existence of an epidemic threshold 
Ac = {k)~^ Uni, such that p = if A < Ac, while p ~ (A — Ac) if A > Ac. In 
this context, it is easy to recognize that the SIS model is a generalization of 
the contact process model, widely studied as the paradigmatic example of 
an absorbing-state phase transition to a unique absorbing state 34 . 

3.2 Uncorrelated complex networks 

For general complex networks, in which large degree fluctuations and corre- 
lations might be allowed, we must relax the homogeneous hypothesis made 
in writing Eq. H13|) and work instead with the relative density Pkit) of in- 
fected vertices with given degree k; i.e. the probability that a vertex with 
k edges is infected. Following Refs. ^3^5) ^^^^ equation for pk{t) can 
be written as 



^The homogeneous mixing hypothesis S3} states that the force of the infection (the 
per capita rate of acquisition of the disease by the susceptible individuals) is proportional 
to the density of infected individuals. 



dp{t) 



p{t) + X{k)p{t) [l-p{t)]. 



(13) 



dt 




(14) 
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In this case, the creation term is proportional to the spreading rate A, the 
density of healthy sites 1 — Pk{t), the degree k, and the variable 6fc(i), that 
stands for the probability that an edge emanating from a vertex of degree k 
points to an infected site. In the case of an uncorrelated random network, 
considered in Refs. |in| lllj . the probability that a edge points to a vertex 
with k connections is equal to kP{k)/ (k) 11 . This yields a 0/^ = Q^'^ 
independent of k that reads as 

^ ' k' 

Substituting the expression ((T^ into Eq. (dJ), it is possible to find the 
steady state solution where Q^'^ is now a function of A alone lllj by the 
following self-consistent equation: 

e- = -L y kp{k) . (16) 

{k) ^ 1 + A/ce-^^ ^ ' 

A non-zero stationary prevalence {pk 7^ 0) is obtained when the r.h.s. and 
the l.h.s. of Eq. (|16|) . expressed as a function of 0'^'^, cross in the interval 
< Q^"^ < 1, allowing a nontrivial solution. It is easy to realize that this 
corresponds to the inequality 



> 1 (17) 

enc=o 



being satisfied. The value of A yielding the equality in Eq. ((T7|) defines 
the critical epidemic threshold Ac, that is given for uncorrelated random 
networks by 

Ar = j^- (18) 

In uncorrelated and infinite SF networks with 7 < 3, we therefore have 
^/c^) = 00, and correspondingly A°"^ = 0. This is a very relevant result, 
signalling that the high heterogeneity of SF networks makes them extremely 
weak with respect to infections. These results have several implications in 
human and computer virus epidemiology |3H]- 



3.3 Correlated complex networks 

For a general network in which the degrees of the vertices are correlated, the 
above formalism is not correct, since we are not considering the effect of the 
degree k into the expression for G^,. This effect can be taken into account, 
however, for Markovian networks, whose correlations are completely defined 
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by the conditional probability P{k' \k). In this case, it is easy to realize that 
the correct factor 0^ can be written as 

Qk{t) = J2P{k'\k)pk,{t), (19) 

k' 

that is, the probability that an edge in a vertex of degree k is pointing to 
an infected vertex is proportional to the probability that any edge points 
to a vertex with degree k' , times the probability that this vertex is in- 
fected, Pk'it), averaged over all the vertices connected to the original vertex. 
Eqs. (|14|) and (|19j) define together the mean- field equation describing the 
SIS model on Markovian complex networks, 

= _p^^t) + Xk[l- pfc(t)] ^ P{k' I k)pk,{t). (20) 

k' 

It must be stressed that this equation is valid only if the network has no 
structure; i.e. the only relevant variable is the degree k. This implies that 
all vertices within a given degree class are statistically equivalent. This is 
not the case, for instance, of regular lattices or structured networks [261 127j 
in which a spatial ordering is constraining the connectivity among vertices. 
The exact solution of Eq. (|2()j) can be difficult to find, depending on the 
particular form of P{k' \ k). However, it is possible to extract the value of the 
epidemic threshold by analyzing the stability of the steady-state solutions. 
Of course, the healthy state pk = is one solution. For small pk, we can 
linearize Eq. pO|l. getting 

^^^L,k' Pk'it). (21) 
k' 

In the previous equation we have defined the Jacobian matrix L = {-L/tfc'} 

by 

Lkk' = -Skk' + XkP{k'\k), (22) 

where 6kk' is the Kronecker delta symbol. The solution /j^ = will be un- 
stable if there exists at least one positive eigenvalue of the Jacobian matrix 
L. Let us consider the connectivity matrix C, defined by C^k' = kP{k' \ k). 
Using the symmetry condition Eq. it is easy to check that if Vk is an 
eigenvector of C, with eigenvalue A, then P{k)vk is an eigenvector of the 
transposed matrix with the same eigenvalue. From here it follows im- 
mediately that all the eigenvalues of C are real. Let Am be the largest 
eigenvalue of C. Then, the origin will be unstable whenever — 1 + AA^ > 0, 
which defines the epidemic threshold 

Ac = T^, (23) 
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above which the solution = is unstable, and another non-zero solution 
takes over as the actual steady-state — the endemic state. 

It is instructive to see how this general formalism recovers previous re- 
sults jlUl lllj . implicitly obtained for random uncorrelated networks. For 
any random network, in which there are no correlations among the de- 
grees of the vertices, we have that the connectivity matrix is given by 
^kk' ~ kP{k'/k) = kk'P{k')/ {k), since the probability that an edge points 
to a vertex of connectivity k' is proportional to k'P{k'). It is easy to check 
that the matrix {C^/).} has a unique eigenvalue = (^k'^'j / {k), correspond- 
ing to the eigenvector v^^ = k, from where we recover the now established 
result Eq. (|TH1) . 

3.4 Correlated Scale-Free networks 

The absence of an epidemic threshold in SF uncorrelated networks is an 
extremely important question that prompts to a possible weakness of many 
real-world networks. Particularly important is for the case of digital viruses 
spreading on the Internet jlUl |1] and sexually transmitted diseases diffusing 
on the web of sexual contacts [HE]- Both these networks show, in fact, SF 
properties that would imply the possibility of major epidemic outbreaks even 
for infections with a very low transmission rate. Immunization policies as 
well must be radically changed in the case that a network has a null epidemic 
threshold [SH EHl EH] • 

In view of the relevance of this framework, it is extremely important to 
study to which extent the presence of correlations are altering these results. 
The main question is therefore which conditions on the degree correlations 
of SF networks preserve the lack of a critical threshold. In the case of cor- 
related networks, we have shown that the epidemic threshold is the inverse 
of the largest eigenvalue of the connectivity matrix C. The absence of an 
epidemic threshold thus corresponds to a divergence of the largest eigen- 
value of the connectivity matrix C in the limit of an infinite network size 
N ^ oo. In order to provide some general statement on the conditions for 
such a divergence we can make use of the Frobenius theorem for non-negative 
irreducible matrices |40j . This theorem states the existence of the largest 
eigenvalue of any non-negative irreducible matrix, eigenvalue which is sim- 
ple, real, positive, and has a positive eigenvector. In our case the matrix 
of interest is the connectivity matrix that is non-negative and irreducible. 
The irreducible property of the connectivity matrix is a simple consequence 
of the fact that all the degree classes in the network are accessible. That 
is, starting from the degree class k it is always possible to find a path of 
edges that connects this class to any other class k' of the network. If this 
is not the case it means that the network is built up of disconnected irre- 
ducible subnetworks and, therefore, we can apply the same line of reasoning 
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to each subnetwork^ . From the Frobenius theorem 001 it can be proved that 
the maximum eigenvalue, A^, of any non-negative irreducible matrix, Ai^^/, 
satisfies the inequality 

Am > inin ^ Akk'i^ik'), (24) 

where {tp{k)} is any positive vector. In particular, by setting A = and 
■0(A:) = k we obtain the inequality 

> min^^fc'^P(^|A:)P(A;'K). (25) 

k' e 

This inequality relates the lower bound of the largest eigenvalue A^ to 
the degree correlation function and allows to find a sufficient condition for 
the absence of the epidemic threshold. By noting that k'P{k' \ i) = 
knn{^i kc), we obtain the inequality 

Al > mmy^iP{i\k)knn{i,kc). (26) 

k ^ — ' 

e 

The r.h.s. of this equation is a sum of positive terms, and by recalling 
the divergence of the ANND with kc, we readily obtain that A^ > oo for 
all networks with diverging (A;^) both in the disassortative or assortative 
cases^. The divergence of A^ implies on its turn that the SIS epidemic 
threshold vanishes, in the thermodynamic limit, in all SF networks with 
assortative or disassortative mixing if the degree distribution has a diverging 
second moment; i.e. a SF degree distribution with exponent 2 < 7 < 3 is 
a sufficient condition for the absence of an epidemic threshold for the SIS 
model in unstructured networks with arbitrary two-point degree correlation 
function. 

In physical terms, the absence of the epidemic threshold is related to 
the divergence of (knn) in SF networks. In homogeneous networks, where 
{knn)j^ — {k), the epidemic spreading properties can be related to the av- 
erage degree. In SF networks, however, the focus shifts to the possibility of 
infecting a large number of individuals in a finite number of contagions. The 
fact that an infected vertex has a very low degree is not very important if a 
hub of the network that provides connectivity to a large number of vertices 
is a few hops away. The infection can, in this case, very easily access a very 
large number of individuals in a short time. In SF networks is the ANND 

^Notice that being irreducible is not equivalent to being fully connected at the vertex 
to vertex level, but at the class to class level. 

^One may argue that, since we are calculating a minimum for k, if the transition 
probability P{£\ko) is zero at some point ko, this minimum is zero. In this case it is 
possible to show that repeating the same argument with C'^ instead of provides us an 
inequality that avoids this problem, Ref. |21|. 
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that takes into account more properly the level of degree fluctuations and 
thus rules the epidemic spreading dynamics. 

It is worth stressing that the divergence of (knn) ^ is ensured by the 
degree detailed balance condition alone, and it is a very general result hold- 
ing for all SF networks with 2 < 7 < 3. On the contrary, the SF behavior 
with 2 < 7 < 3 is a sufficient condition for the lack of epidemic threshold 
only in networks with general two-point degree correlations and in absence 
of higher-order correlations. The reason is that the relation between the 
epidemic threshold and the maximum eigenvalue of the connectivity ma- 
trix only holds for these classes of networks. Higher order correlations, or 
the presence of an underlying metric in the network [^, can modify the 
rate equation at the basis of the SIS model and may invalidate the present 
discussion. 

4 The SIR model 

The susceptible-infected-removed model (SIR) (SB] represents the other paradig- 
matic example of epidemic dynamics. Unlike in the SIS model, in this case 
infected individuals fall, after some random time, into a removed state where 
they cannot neither become infected again nor infect other individuals. This 
model tries to mimic real epidemics where individuals, after being infected, 
acquire permanent immunity or, in the worst case, die. 

The SIR model is defined as follows. Individuals can only exist in three 
different states, namely, susceptible, infected, or removed. Susceptible indi- 
viduals become infected with probability A if at least one of their neighbors 
is infected. On the other hand, infected individuals spontaneously fall in the 
removed state with probability /x, which without lack of generality we set 
equal to unity. The main difference between both models is that whereas 
in the SIS, for A > Ac, the epidemics reaches a steady state, in the SIR the 
epidemics always dies and reaches eventually a state with zero density of 
infected individuals. The epidemic prevalence is thus defined in this case as 
the total number of infected individuals in the whole epidemic process. 

4.1 Uncorrelated homogeneous networks 

In a homogeneous system, the SIR model can be described in terms of the 
densities of susceptible, infected, and removed individuals, S'(t), and 
R{t), respectively, as a function of time. These three quantities are linked 
through the normalization condition 

S{t)+p{t) + R{t) = l, (27) 
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and they obey the following system of differential equations: 

f - 

^ = -p + \{k)pS, (28) 
dR 

These equations can be interpreted as follows: infected individuals decay 
into the removed class at a unity rate, while susceptible individuals become 
infected at a rate proportional to both the densities of infected and suscep- 
tible individuals. Here, A is the microscopic spreading (infection) rate, and 
ik) is the number of contacts per unit time that is supposed to be constant 
for the whole population. In writing this last term of the equations we have 
assumed again, as in the case of the SIS model, the homogeneous mixing 
hypothesis ,23_, 

The most significant prediction of this model is the presence of a nonzero 
epidemic threshold Ac ^^l^. If the value of A is above Ac, A > Ac, the 
disease spreads and infects a finite fraction of the population. On the other 
hand, when A is below the threshold, A < Ac, the total number of infected 
individuals (the epidemic prevalence), i?oo = linit^oo -R(i)i is infinitesimally 
small in the limit of very large populations. In order to see this point, let 
us consider the set of equations (|28|) . Integrating the equation for S'(i) with 
the initial conditions i?(0) = and 5'(0) ~ 1 (i.e., assuming /9(0) ~ 0, a very 
small initial concentration of infected individuals) , we obtain 

5(t) = e-^<'=>^(*). (29) 

Combining this result with the normalization condition (|27|) . we observe 
that the total number of infected individuals Roo fulfills the following self- 
consistent equation: 

R^ = \-e-^'^^)^^. (30) 

While -Roo = is always a solution of this equation, in order to have a 
nonzero solution the following condition must be fulfilled: 

d 



dRr 



1 - e 



> 1. (31) 

floo=0 



This condition is equivalent to the constraint A > Ac, where the epidemic 
threshold Ac takes the value Ac = {k)^^. Performing a Taylor expansion 
at A = Ac it is then possible to obtain the epidemic prevalence behavior 
Roo ~ (A — Ac) (valid above the epidemic threshold). From the point of view 
of the physics of non-equilibrium phase transition, it is easy to recognize 
that the SIR model is a generalization of the dynamical percolation model, 
that has been extensively studied in the context of absorbing-state phase 
transitions |34j . 
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4.2 Uncorrelated complex networks 

In order to take into account the heterogeneity induced by the presence of 
vertices with different degree, we consider the time evolution of the magni- 
tudes Pkit), Sk{t), and Rk{t)i which are the density of infected, susceptible, 
and removed vertices of degree k at time t, respectively El- These 
variables are connected by means of the normalization condition 

Pk{t) + Sk{t)+Rk{t) = l. (32) 

Global quantities such as the epidemic prevalence can be expressed as an 
average over the various degree classes; for example, we define the total 
number of removed individuals at time t by R{t) = X^fc -^(^)-^fc(^)' ^"^^ 
the prevalence as i?oo = linif^oo -R(i)- At the mean-field level, for random 
uncorrelated networks, these densities satisfy the following set of coupled 
differential equations: 



dpkjt) 

dt 
dSkjt) 

dt 
dRkjt) 

dt 



-pk{t) + XkSkm''%t), (33) 
-AA:5fc(t)e'^^(t), (34) 
Pkit). (35) 



The factor Q^^ (t) represents the probability that any given edge points to an 
infected vertex and is capable of transmitting the disease. This quantity can 
be computed in a self-consistent way: The probability that an edge points 
to an infected vertex with degree k' is proportional to k'P{k'). However, 
since the infected vertex under consideration received the disease through 
a particular edge that cannot be used for transmission anymore (since it 
points back to a previously infected individual) the correct probability must 
consider one less edge. Therefore, 

(^E(^-W^)/"^(^)- (36) 

The equations (|33]) , (|3l|) , (|35]) , and (|36|) , combined with the initial conditions 
-Rfc(O) = 0, /Ofc(O) = p^, and S'fc(O) = 1 — p^, completely define the SIR model 
on any uncorrelated complex network with degree distribution P{k). We 
will consider in particular the case of a homogeneous initial distribution of 
infected individuals, p^ = p^. In this case, in the limit p^ 0, we can 
substitute Pfc(O) ~ and S'fc(O) ~ 1. Under this approximation, Eqs. 
and can be directly integrated, yielding 

/•CO 

5,(t) = e-^'=^W, Rk{t)= Pk{r)dt, (37) 
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where we have defined the auxihary function 

m = 1^ G-(r)dr = ^ ^^(fc - l)P{k)R,{t). 



(38) 



In order to get a closed relation for the total density of infected individu- 
als, it results more convenient to focus on the time evolution of the averaged 
magnitude (j){t). To this purpose, let us compute its time derivative: 



d(t){t) 
dt 



(39) 



where we have introduced the time dependence of 5"^ (t) obtained in Eq. 1)3 7|1 . 
Once solved Eq. (jSHJ, 

we Ccin obtain the total epidemic prevalence Roo ^ 
function of (poo = lini^^oo 0(0- Since Rk{oo) = 1 — Sk{oo)^ we have 



(40) 



For a general P{k) distribution, Eq. (|39j) cannot be generally solved in a 
closed form. However, we can still get useful information on the infinite time 
limit; i.e. at the end of the epidemics. Since we have that pk{oo) = 0, and 
consequently limt^oo d</)(t)/dt = 0, we obtain from Eq. (jHU)) the following 
self-consistent equation for (poo- 



1 



(l>-(l>E(*=-i)^'w.- 



(41) 



The value (poo 
the condition 

d 



is always a solution. In order to have a non-zero solution, 



1 



{k) {k) 



l)P{k)e 



-Xk(j>a 



> 1 



must be fulfilled. This relation implies 

±^k{k-l)P{k)>l, 

which defines the epidemic threshold 



A" 



(42) 



(43) 



(44) 



below which the epidemic prevalence is null, and above which it attains a 
finite value. It is interesting to notice that this is precisely the same value 
found for the percolation threshold in generalized networks jSl IH]. This 
is hardly surprising since, as it is well known the SIR model can be 
mapped to a bond percolation process. 
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4.3 Correlated complex networks 

In order to work out the SIR model in Markovian networks it is easier to 
consider the rate equations for the quantities N^{t) and N^{t), defined as 
the number of infected and removed individuals of degree k, present at time 
t, respectively. From this two quantities we can easily recover the densities 
Pk{t) and Rk{t) as 

Pk{t) = ^, Rk{t) = ^^, Sk{t) = l-pk{t)-Rk{t), (45) 

where Nf. = NP{k) is the number of vertices with degree k. The rate 
equations for N^{t) and N^{t) are then given by 



dt 

dN^{t) 



-Nl{t) + XSkit)Tkit), (46) 
Nlit), (47) 



dt 

where we have defined the function 

k' 

In this case, the creation of new infected individuals — the second term in 
the r.h.s. of Eq. H46|) — is proportional to the number of infected individuals 
of degree k' , N^,{t), the probability that a vertex of degree k is susceptible, 
Sk{t), and the average number of edges pointing from these infected vertices 
to vertices of degree k, {k' — l)P{k \ k'), all averaged for all the vertices 
of degree k'. This last term takes into account that one of the edges is not 
available for transmitting the infection, since it was used to infect the vertex 
considered. Dividing Eq. H46|) by N/. and making use of the detailed balance 
condition Eq. ^ we find the rate equations for the relative densities as 



dt 
dpk{t) 



-pk{t) + \kSk{t)ek{t), (50) 



dt 

where, for a Markovian network, the factor 0fc(t) takes the form 

Qk{t) = Y.Pk'{t)^^P{k'\k) (51) 

k' 

Again, it must be stressed that no structure is allowed in the network in 
order for this equations to represent a valid formulation of the SIR model. 

In order to extract information about the epidemic threshold, we proceed 
similarly to the SIS model, performing a linear stability analysis. For time 
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i — > 0, that corresponds to small pk, Rk — and Sk — 1, Eqs. (j5fl]) and (|5T|) 
can be written as 

- ^Lkk'Pk'{t), (52) 

k' 

where the Jacobian matrix L = {Lfcfc'} can be written as 

Lkk' = -Skk' + X ^^'''~^^ Pik'\k). (53) 

In order to infect a finite fraction of individuals, we need the solution pk = 
to be unstable, which happens if there is at least one positive eigenvalue of 
the Jacobian matrix L. Defining the matrix C = {Ckk'}, with elements 

Ckk^ = ^^^^Pik'\k), (54) 

we know from the Frobenius theorem, since it is positive and provided that 
it is reducible at the deg ree class level, that it has a largest eigenvalue 
that is real and positive. Thus, the solution = of Eq. is stable 
whenever — 1 + AA^ < 0. This relation defines the epidemic threshold for 
the SIR model in Markovian networks 

A = J-. (55) 

Am 

In the case of a random uncorrelated network, we have that C^^, = k{k' — 
l)P{k')/ (k). It can be easily seen that this matrix has a unique eigenvalue 
A-m — {k'^) I {k) — ^^ corresponding to the eigenvalue v^'^ = k, thus recovering 
the previous result Eq. (|44|) obtained for this kind of networks. 

The relation between the SIR model and percolation in correlated com- 
plex networks can be closed by noticing that the relevant parameter in this 
last problem is the largest eigenvalue of the matrix C^'^^'^ = {C'^'^f^}, with 
elements C^'j^^ = [k' — l)P{k' \k), as pointed out in Ref. |l2j. It is easy 
to check that if v^^^^ is an eigenvector of C^'^'''^ with eigenvalue A, then 
Vk = k v^^'^ is an eigenvector of C with the same eigenvalue. Then, the 
eigenvalues of C^'^'^'^ and C coincide, yielding in this way the same descrip- 
tion and an identical critical point. 



4.4 Correlated Scale-Free networks 

The discussion of the absence of epidemic threshold of the SIS in SF networks 
with any sort of degree correlations can be easily extended to the SIR model, 
taking again advantage of the Frobenius theorem. In this case, in the general 
inequality given by Eq. ()24() . we set A = and il'ik) = fc, obtaining 

A^ > min ^ Y.^i - l)P{l I k){k' - l)P{k' 1 1). (56) 
k' i 



18 



Given that J2k'i^' I ^) ~ ^nn(^, kc) — l, the previous inequahty reads 

Al > mill V(£ - 1)P{£ I k) [knn{i, kc) - 1] . (57) 

e 

As in the case of the SIS model, the divergence of the ANND with kc in the 
thermodynamic hmit, ensures the divergence of the eigenvalue Am- There- 
fore, a SF degree distribution with diverging second moment is a sufficient 
condition for the absence of an epidemic threshold also for the SIR model if 
the minimum degree of the network is kmin ^ 2. The only instance in which 
we can have an infinite knnii,kc) with a finite eigenvalue Am is when the 
divergence of knn is accumulated in the degree k = 1 and results canceled 
by the term £ — 1 in Eq. (|57j) . This situation happens when the SF behav- 
ior of the degree distribution is just due to vertices with a single edge that 
form star-like structures by connecting on a few central vertices. Explicit 
examples of this situation are provided in Refs. [271 142j . 

5 Conclusions 

In this paper we have reviewed the analytical treatment of the epidemic SIS 
and SIR models in complex networks at different levels of approximation, 
corresponding to the different levels in which degree correlations can be 
taken into account. At the zero-th level, in which all the vertices are assumed 
to have the same degree (homogeneous networks), we observe the presence 
of an epidemic threshold, separating an active or endemic phase from an 
inactive or healthy phase, that is inversely proportional to the average degree 
{k). At this level of approximation, both models render the same result, 
thus showing a high degree of universality. At the first order approximation 
level, in which vertices are allowed to have a different degree, drawn from 
a specified degree distribution P{k), but are otherwise random, we obtain 
epidemic thresholds that are inversely proportional to the degree fiuctuations 
{k"^). The remarkable fact about this result is that the epidemic threshold 
vanishes for SF networks with characteristic exponent 2 < 7 < 3 in the limit 
of an infinitely large network. Finally, in the second order approximation 
level, in which degree correlations are explicitly controlled by the conditional 
probability P{k' \ k) that a vertex of degree k is connected to a vertex of 
degree k', our analysis yields that the epidemic threshold in the SIS and SIR 
models is inversely proportional to the largest eigenvalue of the connectivity 
matrices C^k' = kP(k'\k) and C^k' = k{k' — l)P{k' \ k)/k' , respectively. 
In the case of the SIR model we recover the mapping with percolation at 
the level of correlations exclusively among nearest neighbor vertices. The 
analysis of the divergence of the average nearest neighbors degree knnik, kc) 
with the degree cut-off kc allows us to establish the general result that any SF 
degree distribution with diverging second moment is a sufficient condition for 
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the vanishing of the epidemic threshold in the SIS model. The same sufficient 
condition holds in the SIR model with kmin ^ 2. The SIR model with /cmm. = 
1 always shows the absence of an epidemic threshold with the exception of 
the peculiar case in which the divergence of the average nearest neighbor 
degree is accumulated only on the nodes of minimum degree. These results 
have extremely important consequences, since they imply that correlations 
are not able to stop an epidemic outbreak in SF networks, in opposition 
to previous claims, and indicates that a reduction of epidemic incidence 
can only be obtained by means of carefully crafted immunization strategies 
[S3 OHl inn] ) or trivially through finite size effects jl3] • 
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